Emergent Symmetry and Dimensional Reduction at a Quantum Critical Point 
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We show that the spatial dimensionahty of the quantum critical point associated with Bose-Einstein con- 
densation at T = is reduced when the underlying lattice comprises a set of layers coupled by a frustrating 
interaction. For this purpose, we use an heuristic mean field approach that is complemented and justified by a 
more rigorous renormalization group analysis. Due to the presence of an emergent symmetry, i.e. a symmetry 
of the ground state that is absent in the underlying Hamiltonian, a three-dimensional interacting Bose system 
undergoes a chemical potential tuned quantum phase transition that is strictly two dimensional. Our theoretical 
predictions for the critical temperature as a function of the chemical potential correspond very well with recent 
measurements in BaCuSi206. 
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INTRODUCTION 

The universal properties that appear in the proximity of a 
critical point are determined by a few relevant properties. The 
spatial dimensionality, d , is one of them [1]. This is evident 
from the fact that, in general, the critical exponents depend 
on d. Correspondingly, for strongly anisotropic systems of 
weakly coupled chains or planes, critical behavior character- 
istic for d — \ ox d — 2 respectively can be observed beyond 
a certain distance from the critical point. The critical behavior 
crosses over to three-dimensional only in the close vicinity of 
the critical point of such anisotropic systems. 

In contrast to this conventional dimensional crossover, the 
spatial dimensionality can be effectively reduced under certain 
conditions as the system approaches the critical point. This 
phenomenon of dimensional reduction is closely related to the 
notion of " emergent sliding symmetries"" [3]. Those are phys- 
ical systems for which new symmetry transformations appear 
at low energies (in some cases only at T = 0). In other words, 
the low energy spectrum of the system Hamiltonian is invari- 
ant under these symmetries but the whole spectrum is not [2]. 
We call these transformations "emergent symmetries" because 
they only appear at low energies. By "sliding symmetry" [4] 
we mean symmetry transformations that only change a subset 
of the degrees of freedom which occupy a region of dimen- 
sion lower than d. For instance, if our system is a 3D quantum 
magnet and it is invariant under a spin rotation restricted to a 
given layer, such operation is a "sliding symmetry". 

A simple example of an emergent sliding symmetry is pro- 
vided by classical spins on a body centered tetragonal (BCT) 
lattice with antiferromagnetic XY exchange interactions. If 
the inter-layer exchange interaction, , is smaller than the 
intra-layer one, J||, the energy is minimized when the spins 
are antiferromagnetically aligned on each layer. Since the 
staggered magnetization of each layer can point in any arbi- 
trary direction, the ground state manifold is highly degenerate. 
In this case, an arbitrary spin rotation along the z-axis which 
acts only on the spins of a given layer is a sliding emergent 
symmetry. It is a symmetry because it does not change the 



ground state energy. It is "emergent" because it only exists 
at T = 0: the energy does not remain invariant if we ap- 
ply the same transformation to an excited state. In particular, 
this symmetry is the manifestation of a simple physical prop- 
erty: the order parameters (staggered magnetization) of differ- 
ent layers are decoupled at zero temperature. Consequently, 
in spite of the 3D nature of the system, the antiferromagnetic 
ordering is 2D at T = . This is a simple example of di- 
mensional reduction that results from two key ingredients: the 
classical nature of the degrees of freedom and the frustrated 
nature of the interactions. 
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FIG. 1: Phase diagram obtained by the renormalization group ap- 
proach. 



It is natural to ask if the phenomenon of dimensional reduc- 
tion also exists in quantum systems. In most of the cases, the 
emergent sliding symmetry is removed by zero point fluctua- 
tions. For instance, if we consider now the quantum version 
of the XY model on a BCT lattice with J|| > J^, the ground 
state is no longer invariant under spin rotations along the z- 
axis of all the spins on a given layer. Therefore, this operation 
is an emergent symmetry only in the classical limit. Zero point 
fluctuations remove this symmetry by inducing a finite cou- 
pling between the staggered magnetization on different layers 
[6, 7]. This is a particular example of the phenomenon known 
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as "order from disorder [8]. 

However, zero point fluctuations not always restore the di- 
mensionality by removing the emergent sliding symmetries. 
Such symmetries can appear at special points of the quantum 
phase diagram and lead to dimensional reduction. The main 
characteristic of these special points is that the ground state 
becomes "classical", i.e., it is a direct product of eigenstates 
of a local physical operator For instance, the fully polarized 
ferromagnet is a "classical" state for any spin S. A simple ex- 
ample of dimensional reduction in a quantum system is given 
in Ref.[9] for a Klein model of S=l/2 spins on square lattice. 
In that case, the dimensional reduction from d = 2 to d = 1 
occurs at a first order quantum phase transition point. An im- 
mediate physical consequence of this dimensional reduction is 
the emergence of fractional excitations characteristic of one- 
dimensional systems. 

We have shown recently that the phenomenon of dimen- 
sional reduction can also occur at a quantum critical point 
(second order quantum phase transition) [10]. For this pur- 
pose, we considered the quantum XY magnet of our previous 
example but in the presence of an external magnetic field, H, 
along the Z-direction. The ground state is antiferromgnetic 
for H ~ while the Zeeman term dominates at high fields 
leading to a fully polarized ground state. The antiferromag- 
netically ordered XY component decreases continuously as a 
function of H and vanishes at the critical field He- The spin 
system becomes fully polarized for > He- The correspond- 
ing quantum phase transition is denoted as Bose-Einstein con- 
densation (BEC) because the order is suppressed by suppress- 
ing the amplitude of the order parameter or staggered mag- 
netization. In contrast, the thermodynamic phase transition is 
denoted as XY because the order is suppressed by phase fluc- 
tuations. As we will see below, this difference is crucial for 
the phenomenon of dimensional reduction. The BEC-QCP of 
the system under consideration has a peculiar property: the 
disordered state for H > He is "classical" because it is a di- 
rect product of eigenstates of 5f (z-component of the spin 
operator on a given site i). In other words, the zero point 
phase fluctuations that restore the 3D ordering at iJ = are 
no longer present for H > He simply because the XY spin 
component has been suppressed completely. Since the tran- 
sition is continuous, the 3D coupling induced by these phase 
fluctuations must vanish continuously when H approaches He 
from the ordered side. For this reason, dimensional reduction 
occurs right at the critical point. 

The specific motivation for the theory presented in this pa- 
per is the unusual temperature dependence of the transition 
temperature as function of magnetic field in frustrated mag- 
net BaCuSi206[12]. We describe this system by a Heisen- 
berg Hamiltonian of S* = ^ spins forming dimers on a body- 
centered tetragonal lattice, closely approximating the case of 
BaCuSi206 [13, 14]. The dominant Heisenberg interaction, 
J J2i ^ii ■ ^i27 is between spins on the same dimer i. Since 
there are two low energy states in an appUed magnetic field, 
the singlet and the sf^ + = 1 triplet, we can describe 
the low energy sector either using hard-core bosons, or, in 



terms of the above mentioned XY model. In case of the hard 
core boson description, the triplet state corresponds to an ef- 
fective site i occupied by a boson while the singlet state is 
mapped into the empty site [15, 16]. The number of bosons 
(number of triplets) equals the magnetization along the z- 
axis. The chemical potential ^ = gfisiH — Hd) is deter- 
mined by the applied magnetic field, H, and the critical field 
gfisHci ~ J—2J' (where g is the gyromagnetic factor, /is is 
the Bohr magneton and J' is the inter-dimer exchange inter- 
action). The hoppings ty = J' and t± = (J^ is the frus- 
trated inter-layer exchange interaction) are determined by the 
inter-dimer exchange interactions between spins. Recently, it 
was shown by Rosch and Vojta [17, 18] that the inclusion of 
the two higher triplet modes generates a small coherent sec- 
ond neighbor hopping of low energy triplets between layers 
t*j_2 — J^/J^. This interesting effect restores the d = 3 
character of the spin problem due to the fact that the para- 
magnetic ground state for H < Hd is not purely classical. 
Although it can be described as classical state (direct prod- 
uct of singlets on each dimer) to a very good approximation, 
there are small zero-point phase fluctuations that result from 
virtual process to the higher triplet states (creation and anni- 
hilation of triplet pairs with zero net magnetic moment). It 
was also pointed out in Refs.[17, 18] that the dimensional re- 
duction is still exact at H = Hc2 (saturation field) because 
the state for H > Hc2 is purely classical. For realistic val- 
ues of J = 49.5(1)K and Jj_ < J', J^/J^ < O.lmK in 
BaCuSi206. This implies that the mechanism discussed in our 
paper is still dominant for all experimentally accessible tem- 
peratures T > 30mK. Moreover, the U(l)-symmetry breaking 
terms induced by dipolar interactions will produce a crossover 
to QCP with discrete symmetry at T '--^ lOmK [19] before the 
mechanism of Ref.[ 17] sets in. Finally, the inevitable presence 
of finite non-frustrated couplings in real systems will even- 
tually restore the d = 3 behavior below some characteristic 
temperature Tq < 30mK. 

Despite the above mentioned effects, where lattice distor- 
tions, dipolar couplings or excitations to high energy triplets 
cause a restoration of three dimensional behavior at very low 
temperatures, is it important to stress that the boson model dis- 
cussed in this paper is a nontrivial interacting many-body sys- 
tem where the dimensional reduction at the T = quantum 
critical point is exact. Materials that can be described in terms 
of a chemical potential tuned Bose-Einstein Condensation on 
a frustrating lattice are then candidates for the dimensional re- 
duction as caused by an emergent symmetry in the problem. 
In this sense are the conclusions of our paper are not limited 
to BaCuSi206 alone. 

The main purpose of the present work is to derive the crit- 
ical properties of the field induced BEC-QCP for the XY 
magnet mentioned above. The key finding of our result is 
the detailed phase diagram of Fig.l, where we show the vari- 
ous crossover regimes of a chemical potential tuned BEC on 
a frustrated lattice. This work complements the results pre- 
sented in Ref.[10] by including a renormalization group ap- 
proach (Sec. IV) which provides formal justification for the 
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heuristic mean field approach presented in Ref.[10] and is 
summarized in detail in Section III. The model for the XY- 
magnet on a BCT lattice is introduced in Sec. II. For practical 
reasons, we use the language of hard core bosons which are 
equivalent to 5 = 1/2 spins after a Matsubara-Matsuda trans- 
formation [11]. Our conclusions are presented in Section IV. 



MODEL 

We start from the Hamiltonian of interacting spinless 
bosons on a body centered cubic lattice 

i^B = ^ (Sk - //) a^^Ok + u^n^ni. (1) 



Here rii — ajoi is the local number operator of the bosons 
and 



(2) 



the corresponding creation operator in momentum space. The 
tight binding dispersion for nearest neighbor boson hopping 
on the BCT lattice is 

= £k|| + 2t±7k|| cosfc^c, (3) 

k|| = (kx, ky) refers to the in plane momentum and 

ek|| = i|| (2 + cos k^a + cos kyO) (4) 

is the in-plane dispersion. For convenience we included the 
constant shift 2t\\ in the definition of Eky to ensure that ekj > 
0. The last term in Eq.(3) refers to the inter-plane coupling, 
where the form factor 



7k„ 



cos COS — 

2 



(5) 



describes the k|| -dependence of this coupling in the BCT lat- 
tice. This k 1 1 -dependence is a crucial aspect of our theory. 

For ti_ > and t\\ > t±/2, Bose Einstein condensation 
takes place at Q = (tt/c, 7r/a, A:^). Since 7k|| vanishes for 
k|| ~ (tt/o, tt/o), Eq is independent of k^. The minimum of 
the dispersion is infinitely degenerate as the z-component of 
the wave vector can take any value when the x and y compo- 
nents are equal to ir/a. In case of the ideal Bose gas (u ~ 0) 
this implies for T = that different layers decouple com- 
pletely. Only excitations at finite T with in-plane momen- 
tum away from the condensation point can propagate in the 
z-direction. This behavior changes as soon as boson-boson 
interactions [u > 0) are included. States in the Bose conden- 
sate scatter and create virtual excitations above the conden- 
sate that are allowed to propagate in the z-direction. These 
excitations couple to condensate states in other layers [6]. The 
condensed state of interacting bosons is then truly three di- 
mensional, even at T = 0. This order by disorder argument 
for dimensional restoration due to interactions does not apply 



in case of chemical potential tuned BEC. In this case, the num- 
ber of bosons at T = is strictly zero for /i < 0, i.e. before 
BEC sets in. The absence of particles makes their interaction 
mute and one can approach the QCP arbitrarily closely with- 
out coherently coupling different layers. 

From now on, we will measure the momentum relative to 
the wave vector Qo = (7r/a, tt/o, 0): q = k — Qo, such 
that BEC corresponds to a macroscopic occupation of a state 
with vanishing in-plane momentum, q|| = 0. Since we will 
treat the inter-layer hopping, t±, pertubativly, it is convenient 
to rewrite Hb using real space variables for the direction per- 
pendicular to the planes. 



qii,i x||,!: 
+ t^Y^ 7q||'7y (flqipaqiii + h.c.) , 



(6) 



The indices i,j denote the different layers, with rji,j — 1 for 
nearest neighbor layers while rjij = otherwise. Due to the 
shift of momentum it holds 



sm — sm— . 



7q|| = sm-;:^sm^. (7) 
We note that Hb has a discrete Z2-symmetry [6, 17, 18] 



(8) 



for all i. This is a local Z2-symmetry with respect to the layer 
index. In momentum space the last equation corresponds to 
Qz ^ Qz + tt/c. The in-plane dispersion and the local inter- 
action trivially obey this symmetry. However, the inter-layer 
hopping is only invariant with respect to this transformation 
since 7q|| is odd with respect to either or qy . This discrete 
symmetry is therefore closely connected to the degeneracy of 
the Bose condensate with respect to qz- If we were to include 
an additional inter-layer hopping term between neighboring 
planes with qn -independent hopping ti, 



(9) 



qih'y 



we would break the Z2-symmetry. In addition we would lift 
the degeneracy of the Bose condensate to either kz = tt/c 
or kz = 0, depending on the sign of ti. On the other hand, 
inclusion of a term 



(10) 



that promotes boson hopping between second neighbors 
iVi.i+2 = ^i+2,i = 1 and rjij = otherwise) would lift the 
degeneracy of the bose condensate, but without breaking the 
Z2 -symmetry. We will see below that this leads to impor- 
tant distinction between coherent coupling between nearest 
and next nearest neighbor layers. 
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While the Bose condensed state for /i > and the entire 
regime for T > is three dimensional, the decoupling for 
(/i < 0, T = 0) has dramatic consequences. We show that the 
BEC transition temperature varies as 



where 



Tc (X u In ( — I / In In — 



(11) 



Tc oc /i^/'^ holds instead for an isotropic Bose system in d > 2. 
Despite the fact that different layers are coupled at finite T the 
BEC-transition temperature, Eq.(ll), depends on just like 
the Berezinskii-Kosterlitz-Thouless (BKT) transition temper- 
ature of a two dimensional system [20]. 

The renormalization group (RG) calculation used to ob- 
tain this result (a one-loop RG calculation in analogy to 
Ref.[20, 21] ) shows that the finite temperature transition is a 
classical 3-d XY transition, not a BKT transition. We con- 
clude, therefore, that the T = QCP of chemical poten- 
tial tuned BEC with three dimensional dispersion, Eq.(3), is 
strictly two dimensional. The system then crosses over to be 
three dimensional for /i > or T > 0, where the density 
of bosons becomes finite and boson-boson interactions drive 
the crossover to d = 3. The transition temperature of this 
three dimensional BEC is given by the two-dimensional re- 
sult, Eq.(ll). It is important to stress that the vanishing den- 
sity for (/i < 0, r = 0) implies that these results are not lim- 
ited to weakly interacting bosons[23]. 

MEAN FIELD THEORY 
Phase Boundary 

Here we present a heuristic derivation of Eq.(l 1) based on 
an approach introduced by Popov [25] and further explored by 
Fisher and Hohenberg [20]: infrared divergencies are cut-off 
for momenta g < go — where /i is the chemical po- 

tential. The main results of this approach were already given 
in Ref.[10]. We will show how an effective coupling along the 
z -axis appears when the interaction term of Hb is taken into 
account. For this purpose, we will approach the BEC-QCP 
from the disordered phase. Since we are interested in the case 
of hard-core bosons, we will consider an infinitely large on- 
site repulsive interaction oo. 

While the interaction is local, scattering processes between 
bosons in different layers generate effective non-local interac- 
tions at low energies of the type 



1 



^^int — - 



2 E 

ijfc;,X|| 



v^jki (X||)4,|j4|,jax||fcax,,i- 



(12) 



To leading order in the boson density p, the Fourier transform 
wo(q|| ) = Van (q|| ) of the intra-layer effective interaction re- 
sults from adding the ladder diagrams shown in Fig. 2a [26]: 



1 



i^o(qi|| + q2||) 



1 



47r2 e 



(13) 



t|| (2 — cos qxO, — cos qyo) 



(14) 



due to the shift of the in-plane momentum. The integral 
in Eq.(2 ) diverges logarithmically in two dimensions for 
qi||jq2|| — ^ 0. The effective interaction will be logarithmi- 
cally small in the low density limit. An heuristic way of de- 
riving a consistent mean field theory is to introduce the cut-off 
qo^ VM7i[20, 24, 25]: 



1 



1 



In ^ 



47r2 ec 



(15) 



(a) 



_y n.k,+k3 



n,k 



n,k2-k:f 



n+l,ki+kj n+l,ki-kj 




n+2,k,+kj n+2,ki-kj 




FIG. 2: (a) Ladder diagrams that provide the dominant contribution 
to the intra-layer scattering in the low density regime [26]. (b) and 
(c) leading order diagrams that contribute to the coherent inter-layer 
hoppings tx.i ™d il,2- 

We proceed now to compute the inter-layer interactions 
Vijki that are generated by combining the intra-layer renor- 
malized interaction, v^, with the inter-layer hopping term i^. 
For this purpose, it is convenient to expand the propagator in 
powers of the inter-layer hopping: 



Gi3i.n) = 9{q\\)5i] +Uj{q\\)g'^{q\\) 
I 

where Uj (q|| ) = t_L7q|| Vid and 
9{<1\\) 



(16) 



(17) 



is the intra-layer propagator for long wavelengths (g|| <C 1). 
From now on, we measure in plane momenta in units of 27r/a 
where a is the in-plane lattice constant (i.e. we set a = 1) and 
work in the long wavelength limit g|| ^ 1. The leading order 
inter-layer effective interactions that are relevant for inducing 
coherency along the z -axis are vujj = for \i — j\ 



q™ +qiii 



1 (see Fig. 2b) and 



2 (see Fig. 2c). Analyzing the 
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corresponding ladder diagrams yields: 

dujcPq 



2,2 



(2^) 



lqlp-q9 ( [p-q), 



Tqip-qg iq)9 ip-q)- 



Performing the momentum and frequency integration with 
lower momentum cut off go and setting p yields 



2-l2 



Vl 



V2 



o In — , 



2+4 



9v^t 



■In 



1287riji qo ' 



(18) 



The vq, Vl and V2 processes generate the minimal number 
of terms that have to be included in the low energy effective 
Hamiltonian in order to provide a correct description of the 
critical properties of our bosonic system in the low density 
limit. The expression for the new interaction term in the low 
energy theory is 

I 2 

^int 2 ^ VrnSii^jim ial,, x,, jaxn j • (19) 

X|| ij m—0 

The m = term corresponds to intra-layer scattering ver- 
tex Vq- The other terms with m = 1, 2 describe hopping of 
pairs of bosons from the layer i to the layer j ± m. Now we 
perform the mean field decoupling (we suppress the in-plane 
coordinate X|| for simplicity) 



where p 



With this mean field approximation we 



obtain an effective single particle Hamiltonian with disper- 
sion 



i? * = i?q + 2viKi COS qz + 2v2K2 COS 2qz, (20) 



with 



4,^2 \"q||,i"q||,J+i/ ' 

and effective chemical potential 

/i* = ^ - vop. 
The mean values (al m„ , ,, ) are given by: 



(21) 



(22) 



"q||,i"q||.i+J 



dqz cosjjq^) 

27T AK'^^') - 



1 



It follows (flqii jOqii = 0, a rcsult that is a consequence 
of the local Z2 symmetry of H. This means that ki may only 
become non-zero when the U{1) symmetry gets broken at the 



BEC transition. In contrast, {a]^^^ ja^n ^^2) 1^ invariant under 
the discrete Z2 -symmetry of H. Therefore this mean value is 
finite as long as the concentration of bosons is finite. Although 
the term 2t±jk^^ cos qz cancels at = , it is crucial to keep 
it in order to obtain a finite value for K2- Without this term, we 



.a, 



qii ,n+2/ 



have: i?* = which would imply (a^^ 

0. The system undergoes a Bose-Einstein condensation when 
the effective chemical potential, ji*, becomes equal to zero: 



/i = vop{Tc). 



(23) 



In order to calculate p{Tc), we need to solve the integral 
(21) for K2 at T = Tc. We will assume that Sc^ > 
2v2K2 cos2(jz for any q\\ > qo, and evaluate the expectation 

value (^ajjii n'^c^^^ n+2^ the limit K2 = 0. Below we verify 
that this assumption is justified for small t±_ /t\\ . It follows 



K2 



Tiling 



(24) 



where the logarithmic term contains again the lower momen- 
tum cut off qo ■ Without this lower cut off the mean field 
theory could not be properly defined. This result is con- 
sistent with the above assumption that 2v2H2 is small com- 
qo, since Sq^ ~ typ/ ln/.i/f|j while 

iQ. The last result was obtained using 



> 



pared to equ if g|| 

2V2K.2 ^ (^i/^ll)*^ 

Eq.(18)foru2. 

With K2 7^ for finite T, the effective dispersion i?* 
of Eq.(20) becomes three dimensional. Coherent motion 
of bosons within the planes and between planes is allowed. 
While thermally excited bosons are needed for this coherent 
hopping to emerge, its origin are quantum fluctuations that 
cause the non-local interlayer interaction W2- The quantum 
critical point at (T = 0, yU = 0) is however purely two dimen- 
sional. We have a finite coherent inter-layer coupling at the 
BEC momentum only for finite T or in the bose condensed 
state. This implies that the bose condensate itself is three di- 
mensional and that the universality class of the finite T tran- 
sition is 'dD — XY. However, the amplitude of this coherent 
coupling is very small and the system will be effectively two 
dimensional until it is very close to the transition. The width 
of the regime with three dimensional fluctuations shrinks to 
zero as Tc vanishes. This implies that the magnitude of Tc 
obtained from Eq.(23) is practically the same as the magni- 
tude of the Kostelitz-Thouless temperature Tkt- The thermo- 
dynamic phase transition is however always of second order 
Since the effective coupling V2 induced by order from disorder 
is irrelevant for the quantum critical point (the phase transi- 
tion induced by changing the chemical potential at T = 0), it 
is relevant for the classical phase transition at T = Tc- There- 
fore, the dependence of Tc on 1^1 and p is given by the d = 2 
expressions: 

P 



Tr OC ill 



p OC 



In In p~ 
Tc 
InU/Tc' 



(25) 
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In addition we have the usual two dimensional expressions 
for the density as function of temperature or chemical poten- 
tial: 



p{T = 0,fi) (X (Uln— , 
T 

p{T,^i = 0) oc — In In ill /T, 



(26) 



The appeal of this mean field theory is its physical trans- 
parency and technical simplicity. The introduction of the 
chemical potential as lower cut off is however rather ad hoc 
and it is unclear whether it is justified for the problem at hand. 
In order to avoid these ambiguities we developed a renormal- 
ization group approach that confirms the results of this section 
(see below). 



Bond Ordering 

We will discuss now the bond ordering that accompanies 
the BEC. The Z2 -symmetry (8) is spontaneously broken be- 
low Tc because according to Eq.(21) 



1 

4^' 



-.t 



(27) 



becomes finite for a nonzero BEC order parameter (flg j). 
Moreover, | ki | is maximized when the relative phase between 
(oq j) = Ae"*' and (cq ^^i) = Ae'^'+'^ is or tt meaning that 
the inter-layer coupling favors any of these two relative ori- 
entations below Tc'. oc cos^ {4'i+i ~ [6]- In real 
space, this means that the phase of a given site X|| of the layer 
i, 0X|| ,i, is parallel to the phase of two of its nearest-neighbors 
on layer i + 1 and antiparallel to the phase of the other two. 
Consequently, the four bonds connecting a given site with its 
nearest-neighbors on an adjacent layer become inequivalent 
below Tc, i.e., there is a finite bond order parameter 

In principle, the bond ordering could appear at a critical 
temperature higher than Tc- In that case there would be two 
thermodynamic phase transitions instead of one. We will 
show now that there is only one phase transition, i.e., that the 
bond order parameter becomes continuously nonzero only be- 
low Tc- For this purpose, we introduce: 



(28) 



According to Eq.(20), the transition to the Bose condensed 
state occurs for = 0. Therefore, 5p measures the deviation 
of the chemical potential from its critical value. By computing 
the integral (21) for j = 1 we obtain: 



T 



(1 - y)dy 



In I - e.-iP'^ilv+P^t^+Pi^/H) 
(29) 

where t\ = — vi^i and we have used the heuristic cut-off qq. 
We do not expect any transition for Sp,/ fi ^ 1 because the 
temperature is much smaller than the excitation gap and the 



number of bosons becomes exponentially small. Therefore, 
we will assume Sfi/fi <C 1 that corresponds to the quantum 
critical regime with the temperature (or the chemical poten- 
tial) approaching the BEC-point from the disordered side. If 

t*^/n <C 1, we obtain: 



T 



Kl = 



(1 - y)dy 
that reduces to 



III if32tly + f3Sfi + f3ij/t\\), 
(30) 



^ - 1 _ 
fi 47r/it|| 



(31) 



after expanding the logarithm. Eq.(31) violates the original 
assumption ^ 1 meaning that there is no solution of 

Eq.( 29) for any finite This imphes that the bond ordering 
appears only below Tc- 



RENORMALIZATION GROUP APPROACH 

The mean field approach presented in the previous section 
was supplemented by the introduction of a lower cut-off of 
otherwise infrared divergent terms in the perturbation theory. 
These divergencies result from the fact that the two dimen- 
sional dilute Bose system is a quantum system at the upper 
critical dimension. The natural approach to control these di- 
vergencies is a renormalization group analysis. 

In our renormalization group analysis of the model Eq.(l) 
we start from the action: 



Share = X! / * 



'9J 



(32) 



E 



where 



i •'91,2,3 



Goij (<?) = 9 ^(9)% - t^qxQyn, 



(33) 



with g{q) of Eq.(17). Here q = (q||,cj„) refers to the pla- 
nar momentum q|| ~ {q^, qy) and the bosonic Matsubara fre- 
quency LOn — 2nnT. We use the notation 



T 



q|||<A (27r) 



(34) 



The upper cut off A is determined by a length scale larger than 
the inter-atomic spacing but much smaller than the correla- 
tion length. Thus A ~ 1 with our choice that the in-plane 
lattice constant a = 1. The upper momentum cut off A yields 
an upper energy cut off of order . 

Although the inter-layer hopping tx_ is a marginal pertur- 
bation, it is responsible for the emergence of new, non-local 
interactions, where excited states of one layer propagate into 
another layer and couple to its low energy states as it was al- 
ready demonstrated in the mean field theory of the previous 
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section. We need to include such non-local couplings into th( 
effective action of the renormalization group analysis. Sucl 
non-local interactions might cause, in turn, coherent motioi 
of bosons even for q|| — 0. Thus, we need to further supple 
ment the action and allow for a coherent motion of bosons 
This leads to the effective action: 



J2 / 



4e 



(35 



ijkl "'919293 



where 



(9) =5 -U3 (q|i) 

with inter-layer hopping 



(36 



Thus, we include terms like Ti 2 in Eqs.9 and 10. In particular, 
the hopping t2 between next nearest neighbors is included as 
it is the leading inter-layer boson hopping that does not violate 
the above mentioned Z2-symmetry. The hopping ti between 
neighboring layers is included to explicitly demonstrate that 
it will not contribute to coherent interlayer motion. At the 
beginning of the renormalization group flow S = ^barc and 
it holds that 



and 



h (1 = 0) = t2 (1 = 0) = 0, 



Vijki {I = 0) = uSijSikSii. 



(37) 



(38) 



Before we derive one loop RG equations we discuss the var- 
ious distinct physical regimes indicated in Fig. 3. The renor- 
malization group approach is controlled by the flow variable 
I — log (A/A (Z)). As usual in the regime of small but finite 
temperature, T is a relevant perturbation of the T = QCP 
and a crossover to classical critical behavior occurs for I > Iq 
when the renormalized temperature becomes comparable with 
the upper energy cut off t|| : 



(39) 



Excitations in the system with momentum larger than Ae^'" 
behave just like T = quantum excitations, while those 
with momentum below Ae"'" are classical. As T — ^ 0, the 
crossover variable Iq ^ 00 and, as expected, all degrees of 
freedom are in the quantum regime. 

In addition to this quantum to classical crossover, the sys- 
tem undergoes a dimensional crossover at a scale I* defined 
via 



{I* 



(40) 



depending whether ti (l*) or ^2 (l*) first reaches t\\. In anal- 
ogy to the quantum to classical crossover, it holds that exci- 
tations with momentum larger than Ae^' behave quasi two 



pair scattering 
between layers 



weakly coupled 
layers 



isotropic 
3D behavior 



2D quantum 2D classical 

regime regime 



3D classical 
' regime '=^'>B(^Mt)) 



FIG. 3; Regimes of the renormalization flow. 



dimensional while those with momentum below Ae"' are 
sensitive to a coherent inter-layer coupling. We show that 
lo < I* if the system is close to the critical temperature, i.e. 
the dimensional crossover is driven by the existence of ther- 
mal excitations in the system. However, quantum fluctuations 
are nevertheless crucial for the dimensional crossover, as they 
lead to non-local interactions Vijki that are responsible for the 
dimensional crossover once thermally excited bosons exist. 

As long as the system is in the regime I < I* , the renor- 
malized coherent inter-layer coupling is small. This has im- 
portant implications for the distinction between low and high 
energy degrees of freedom. For I < I* we have to integrate 
out states with Ae~' < |q|| | < A, regardless of the momen- 
tum perpendicular to . Only once the RG flow enters a 
three dimensional regime for I > /* is it sensible to distin- 
guish low and high energy modes with momentum per- 
pendicular to the planes. Then we integrate out states with 
-/ 



A 



< A. 



qji 

We first give the one loop renormaUzation group equations 
fovl <r . It holds 

-j7 = 2// - 2 ^ vumi I Gi„i (q) , 

Irn ■^f' 



dt± 
~dl 
dt 



^ = 0, 



=■ 2tj + 2^ Wij^,„^i+j / Gi„i{q), 



dT 
Hi 



2T, 



(41) 



as well as 



dv 



ijlm 



dl 



stuv "^9 

4:'^VisiuVjtmv Gsv{q)Gtu{q) 



stuv 

where we use the short hand notation: 



/ ... = limZ-^T V / d'^q\i.. 

Jq „ JAe-'<|q|||<A 



(43) 



For I > I* the renormalized inter-layer hopping element 
is comparable to the in-plane hopping and we are finally 
allowed to perform a continuum's theory for the direction per- 
pendicular to the layers as well. Then, the problem is identical 



7 



RENORMALIZATION GROUP APPROACH 



to the one of an isotropic three dimensional Bose system 



Viso {I*) aQ^aQ^aQ3aQi+Q2-Q3 i^) 



where Q = (q|| , qz,0Jn) is a (3 + 1) -dimensional vector that 
includes the momentum Qz- The initial values for this flow are 
determined by the final values for the flow for I < I* . The 
isotropic boson interaction 



C^Vijkl {I*), 
ijkl 



(45) 



corresponds to the Qz = value of the coupling constants 
Vijki- For short ranged couplings Viso it is dominated by 
the largest coupling constant. The additional prefactor c, 
with lattice constant in the z -direction, ensures that viso 
is a three dimensional coupling constant with dimension 
[length]"^ [energy] for d = 3. 



Quantum-Classical Crossover 

We first analyze the quantum to classical crossover at I — 
Iq. It is useful to consider separately the behavior in the 
quantum regime, where the renormalized temperature is small 
compared to the upper energy cut off T (l) < ty and the clas- 
sical regime, where T (1) > We treat both regimes sepa- 
rately and assume T (l) ^ <|| in the former and T (l) S> t|| 
in the latter regime and connect the flow of the various cou- 
pling constants smoothly at Iq. This is essentially the approach 
taken in Ref.[20]. The analysis of Ref.[22] demonstrates that 
the approach used here is fully consistent with results obtained 
using a more careful analysis of the crossover behavior. 

At T = 0, the flow of the chemical potential and of the co- 
herent inter-layer boson hopping are unaffected by the in- 
teraction between bosons. This result is specific for the prob- 
lem of dilute bosons with i^l < since for T = 



Glm (q) = 0, 



(46) 



as a result of the integration over frequency. Physically this is 
due to the fact that the boson number vanishes for T = Q and 
/i < 0. This yields the renormalization group equations; 



diJL{l) 
dl 

dti,2 (0 
dl 



Ml), 
2tl,2 (0 . 



(47) 



At T = 0, interaction do not affect /i, ti and t2. This only 
happens once the system reaches the classical regime I > l^. 
Since ti 2 (0) = 0, it follows that ii 2 (0 = as long as the 
system is in the quantum regime. Quantum fluctuations do 



not induce a coherent hopping between layers. The inter-layer 
hopping (k||) — Tyiji^fc^fcy remains unchanged. Note, that 
the amplitude t± of this inter-layer coupling is unchanged un- 
der renormalization. The flow of the chemical potential is 



= /ie 



2/ 



(48) 



Next we analyze the behavior of the interactions. At T 
holds 



Gsv (q) Gtu (q) = 0, 



(49) 



which vanishes again because of the vanishing Boson density. 
We are left with the analysis of 



dv. 



ijlm 



dl 



IJUV 



Vstlr. 



Gsu (q) Gtv i-q) . (50) 



In Appendix A we analyze this flow equation in the limit 
where the bare inter-layer hopping t± is much smaller than 
the in-plane hopping t||. Up to order (t.±/t\\) we have to an- 
alyze the coupling for bosons in the same layer vq = vuu , 



in neighboring layers vi 
neighbor layers vi = vuj 
I (I > 27rt||/u): 



= Viijj with j — i ± 1 and second 
with j = i ± 2. We obtain at large 



vo (0 
V, (l) 

V2 (0 



7r<||a 
2/ 

Trill 5 / , ^4 



(51) 



It is important to keep in mind that these results were obtained 
with the assumption that initially uo = 0) ~ u is the only 
coupling constant. The inter-layer interactions v\ and V2 re- 
sult from multiple scatterings in distinct layers where virtual 
bosons propagate between layers. Thus, we find that there is 
no coherent coupling between layers in the quantum regime 
i.e t\,2 (0) = 0. On the other hand, we do find that non-local 
interactions, that couple different layers, emerge. This is fully 
consistent with the finding of Refs.[7? ]. 

At finite T, the flow in the quantum regime stops at 



(52) 



For / > Iq thermal, as opposed to quantum fluctuations, come 
into play. The initial values for the subsequent flow are of 
course the final values of the RG flow of the quantum regime: 
[I (Iq) = /ie^'" = t\\ ^ and Vi (Iq) where the Vi (/) are given in 
Eq.(51). 



Dimensional Crossover 

The RG flows for I > Iq continues to be two dimensional, as 
no coherent inter-layer was generated in the quantum regime. 
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As discussed above we will now analyze the flow equations 
as if the problem was purely classical, i.e. we include solely 
the lowest Matsubara frequency in the evaluation of the Feyn- 
man diagrams. In this case temperature only enters the flow 
equations in the combination 



where 



Wi (0 



T{1) 



V, (0 . 



Thus, it is convenient to use m; (1) in what follows. The lead- 
ing order flow equations of Wi (l) are 



dwj (l) 
dl 



2w, (0 ., 



with solution 



(54) 



(55) 



The coupling constants Wi [l) are relevant. This is a conse- 
quence of the fact that the upper critical dimension of the 
classical regime is du. class. = 4 as opposed to du.qu. = 2 
for the zero temperature quantum regime. If we wanted to de- 
termine the critical exponents of the classical phase transition, 
we would have to include higher order terms. As pointed out 
by Millis[22], it is not necessary to include these higher order 
terms if one only wants to determine the value of the transition 
temperature: At low T, the coupling constants Vi{l) decrease 
for large Iq as follows from Eq.(5I). Thus, the initial values 
Wi (Iq) of the classical flow are small. While the interactions 
become relevant for I > Iq corrections to Eq.(54) remain neg- 
ligible unless the flow enters the actual critical regime. How- 
ever, in our case the flow only enters the critical regime after 
the dimensional crossover. Thus, we can, for the moment, 
safely neglect corrections beyond Eq.(54). 

As shown in Appendix B, the RG flow equations for the 
coherent hopping elements and the chemical potential in the 
classical regime are 



dfi 

m 

dti 

~dl 

dh 
dl 



2fi Wo, 

TT 



t2 

2t2 + W2. 

■Kti 



(56) 



It immediately follows that ti {I) = since ti {Iq) = 0. No 
coherent nearest neighbor hopping ii is being generated by 
the mechanism we describe. This is a consequence of the dis- 
cussed Z2-symmetry. A finite value for ti corresponds to a 
broken Z2 symmetry. However, the second neighbor coupling 
t2 flows to a finite value even if its initial value vanishes. If we 
use Wi (l) of Eq.(55) with initial values Vi (Iq) from Eq.(51) it 
follows 



g2(i-io) 



dl Iq 
dl Iq 



(57) 



1 3 

90 = 4(1- 2 (i±A||)'+32 (^^All 



92 



(58) 



(53) The solutions of these differential equations are 



m(0 

t2 (l) 



m('o) 



{I - lo) 



2(l-lo) 92t\\ jl-lp) 



la 



(59) 



In the last equation we already took into account that the initial 
value of the coherent hopping vanishes: t2 (Iq) = 0. The 
dimensional crossover takes place at I* where \t2 {l*)\ — 
which corresponds to 



2(i'-io) 92 jl* - Iq) 
For large I* — Iq this is equivalent to 



1. 



92 log /0/52 

This yields 

wo {n=vo {lo)e^^''-'°'> =vo{lo) 
as well as 



lo 



92 ^ogk/92 



lo 



.92 log (0/32 92 



Inserting Iq ~ ^ log (eo/T) and fi (lo) gives 
5o7r<|| 1 



wo in 



2 52log(ilog(t||/r)/.92) 



(60) 



(61) 



(62) 



(63) 



(64) 



for the value of the coupling constant at the end of the two 
dimensional flow and 



nil*) = i|| 



-2^os{h/T) 



a2log(ilog(t||/T)/g2)T .92 



^-^ (65) 



for the corresponding chemical potential. As pointed out 
above, for I > I* , the RG probes energies sufficiently low to 
be sensitive to the three dimensional character of the system. 
These final values of the combined quantum and classical two 
dimensional flow become the initial value of the three dimen- 
sional flow. Since always holds /* > Iq, \i follows that this 
three dimensional flow is always in the classical regime. 



Flow in the Three-Dimensional Classical Regime 

The final regime of the flow is in the classical three dimen- 
sional regime. The flow equations are the usual ones for an 
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isotropic three dimensional classical bosonic system, i.e. for 
an two component ip'^ or XY model. The condition for the 
critical temperature is that the initial values for the flow of 
this three dimensional classical flow obey: 



(66) 



This ensures that the flow is on the critical surface and the 
system is close to the critical temperature. An alternative way 
to interpret this condition was given in Ref.[22] where it was 
shown that Eq.(66) is equivalent to the Ginzburg criterion for 
the onset of critical fluctuations. Whenever a system is in the 
Ginzburg regime of classical critical fluctuations, it is very 
close to the actual critical temperature. The detailed analysis 
inside this regime is the usual one for a d = 3 classical XY 
model and does not need be reproduced here. We are more 
interested in the value of the transition temperature at low T. 
We use our previous results for the initial values /i (l*) and 
w {I*) of the three dimensional flow to analyze the condition 
Eq.(66) and obtain 



pi = Tc 



2go(f +log (4 log (t||/T,; 



ST, 



log(£o/T;) 
log(^log(t||/r,)) 
log (eo/T,) 



(67) 



Solving this result for Tc with logarithmic accuracy yields the 
transition temperature as function of chemical potential, as 
given in Eq.( 11). The phase diagram that results from our 
RG analysis is represented in Fig.l. 



SUMMARY 

In summary, we have shown that inter-layer frustration re- 
duces the effective dimensionality of a BEC quantum phase 
transition induced by a change of the chemical potential. The 
BEC-QCP exhibits 2D quantum critical fluctuations that dom- 
inate over an extended region of the phase diagram. The phase 
boundary between the disordered and ordered phase extends 
to finite temperatures although the universality class of the 
transition changes from BEC in 2h-2 dimensions at T = to 
3D-XY at finite T. For T > 0, there is a crossover from the 
2D quantum critical to a 2D classical regime as the system 
approaches the phase boundary from disordered side. The di- 
mensional crossover occurs within the classical regime as the 
system gets even closer to the phase boundary (see Figs.l and 
3). 

The BEC ordering is accompanied by bond ordering that 
results from a spontaneous breaking of the Z2-symmetry dis- 
cussed in section 11. Both, the BEC and the bond order pa- 
rameters increase continuously from zero for i^l /ic. A fi- 
nite bond-order parameter induces a finite hopping between 
nearest-neighbor layers that vanishes at the phase boundary 
together with the bond ordering. 

Although according to our results the thermodynamic phase 
transition always belongs to the 3D-XY universality class. 



this transition becomes more quasi-2D like as the system ap- 
proaches the 2D BEC-QCR Besides the consequences that 
were already discussed in the paper, like the peculiar behav- 
ior of Tc(/i) given by Eq.(ll), this observation has implica- 
tions for the T-dependence of any thermodynamic quantity 
for T ^ Tc and /i > iic{T = 0) given the dimensional 
crossover predicted by our RG calculation. 

The dimensional reduction at a QCP that we discussed in 
this paper can be experimentally verified in real quantum mag- 
nets such as BaCuSi206 [10, 12]. For quantum magnets, 
the chemical potential corresponds to a magnetic field ap- 
plied along the symmetry axis while the particle density cor- 
responds ot the magnetization per site. Therefore, the quan- 
tum phase transition discussed in this paper corresponds to 
the suppression of magnetic XY-ordering by the application 
of a magentic field that saturates the moments along the Z- 
direction. Although we discussed the case of a BCT lattice, 
our result can be applied to more general layered structures 
with frustrated inter-layer coupling. 
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APPENDIX 

Flow in the 2D-quantum regime 

In this appendix we derive the result, Eq.(51) for the inter- 
actions Vi in the regime I < Iq prior to the quantum to classical 
crossover, by solving flow equations 

dv^jim _ f Gsu{k)Gtv{-k) . (68) 



dl 

We derive these results by expanding with respect to the ra- 
tio of the hopping elements perpendicular and parallel 
to the layers. Thus we expanded the propagator Gij (k) of 
Eq.(36) in powers of the inter-layer hopping, see Eq.(16 ). In 
perturbation theory in tj_, it always holds that i = j, I ^ m , 
u = V , and s — t. Thus we obtain (to simplify the notation 
we use Vij = vujj): 

St ''1 

Including terms up to order t"^ it follows: 

S st 

-Y,vuVsjJ2vsimtA<^^\ (70) 
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where 



liiTgii) ' 5 (-9)" 



9 

2 (•> 



with 7 (fc) = fcj-fcy and 5 (q) of Eq.(17). Performing the 
frequency and momentum sums yields = (27rt||)-\ 

= ^(87ri||)"\ ^ ^(2)/2 as well as A^^) = 

II 

tytj|x9 (l287r<||)"\ Only terms with j = iilandj = i±2 
are being generated at fourth order in t^. We will then intro- 
duce three different coupling constants, wq = vn, Vi = Vi i±i 
and V2 ~ Vi,i±2- It will turn out to be crucial to include 112 in 
addition to the leading non-local coupling vi . Performing the 
lattice sums yields explicit flow equations for the three cou- 
pling constants. If we now keep in mind that due to the initial 
conditions vi {I = 0) = V2 {I = 0) ~ vertices with vi are at 
least of order t']_ and vertices with 112 are at least of order <^ 
we can restrict the flow equations to fourth order in t±: 



^ = -2t;2A(o)-4z;ot;iA(2), (71) 

with = y4(") + 4B(2) + 2A^'^'>. For large ^ we expect a 
decay of the coupling constants according to Va (l) oc l^^ . 
Thus we assume 



MO 
I ' 



(72) 



and analyze the flow equations for (l). For large enough 
I we can determine the amplitudes of the coupling constants 
from = 0, leading to the algebraic equations 

ho = hlA^°'> + 2hjA^°'> + 4hohiA^^\ (73) 
hi = 2Mi(a(")+4B(2)) +hlA^^\ 

h2 = {2hah2 + hi) A'-^^ +2hahiA'-^^ +hlA'^'^\ 

We can solve this system of equations once again by expand- 
ing with respect to the small parameter 



(74) 



keeping in mind that A^^^ /A'^^^ = |a^. It foflows 



ho = 2TTt»a^ 1 - 2a 



hi = -27rt||a^ (a ~ 4a^) 
/i2 = —2TTt\ia^—a'^. 



(75) 



Flow equations in tlie 2D-classical regime 



As discussed in the main text, in the two dimensional clas- 
sical regime we concentrate of the flow equations of the chem- 
ical potential and coherent hopping elements. We start from 
the general RG equations given in Eq.(41). Using the fact that 
Vijki has only three nonvanishing contributions Vm = vnjj 
with j ^ i ± m and m = (same layer), m = 1, neighbor- 
ing layers and m = 2 (second neighbor layers). Inserting this 
result into Eq.(41) yields: 



^ = 2^^-2vo r Gu{k), 
dl Jk 

2ti+2vi [ Gu+i{k), 

Jk 



dti 
H 

dl 



2^2 + 2^2 / G«»+2(fc). 
Jk 



(76) 



It holds up to second order in t^: 

Gu{k) ~ g{k) + 2tl^{kfg{kf, 
Gii+i{k) ~ ti_^{k)g{kf , 
Gu+2{k) ~ tlj{kfg{kf. 



(77) 



Here we ignored effects due to ii and t2 as those will only be 
of higher order in t± /i || . This enables us to perform the shell 
integration 



Gii (fc) 




G.,+1 (fc) 



Gu+2 (k) 



0, 



Tjl) t\ 
2-KaH\\ tl 



(78) 



where we only included the zeroth's Matsubara frequency in 
the classical regime. The contribution for the nearest neigh- 
bor coupling vanishes since J^, 7 (fc) = 0, an effect caused by 
the Z2 -symmetry of the Hamiltonian. Inserting these results 
into Eq.(76) yields Eq.(56). The solution of the flow equations 
then yields values for the second neighbor hopping small by 
(tj^/^ll)^, justifying our assumption to neglect 1.2 in the right 
hand side of Eq.(77). It is also important to notice that includ- 
ing terms with coherent neighbor hopping t\ in Eq.(77) and 
self consistently solving the RG equation for ii still yields 
ii = on the disordered side of the phase transition. 



Inserting these results into Eq.(72) yields the result Eq.(51 ). 
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